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Abstract 

We study a series array of nonlinear electrical circuit elements that possess negative differential 
resistance and find that heterogeneity in the element properties leads to the presence of multiple 
branches in current- voltage curves and a non-uniform distribution of voltages across the elements. 
An inhomogeneity parameter r max is introduced to characterize the extent to which the individual 
element voltages deviate from one another, and it is found to be strongly dependent on the rate of 
change of applied voltage. Analytical expressions are derived for the dependence of r max on voltage 
ramping rate in the limit of fast ramping and are confirmed by direct numerical simulation. 

PACS numbers: 72.20.Ht, 84.30.-r 
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Negative differential resistance (NDR) - in which an increasing applied voltage causes a 
reduced electrical current flow - occurs in a number of electronic transport systems, for exam- 
ple, tunnel diodes,- semiconductor double barrier quantum well structures,- and molecular 
bridges.- 1 ^ When several such elements are placed in a series array, forming a complex system, 
the NDR of individual elements often leads to striking self-organization effects associated 
with a spatially non-uniform distribution of electric field, despite the fact that the elements 
are almost identical. An important example of this behavior is provided by semiconductor 
superlattices consisting of a series array of quantum wells (e.g., composed of GaAs) sepa- 
rated by potential barriers (e.g., composed of AlAs). In the case of doped superlattices with 
relatively wide barrier layers (i.e., weakly-coupled superlattices), one may observe multiple 
current branches in current-voltage (J — V) characteristics with the number of branches 
approximately equal to the number of periods of the superlattice.-^ These branches are 
associated with a static, non-uniform electric field configuration in which some periods 
of the structure are in a low-field domain and the others in a high-field domain.- 1 ^ Cur- 
rent branching and/or non- uniform electric field distributions have been reported for other 
spatially-periodic systems such as quantum cascade laser structures,- multiple-quantum well 
infrared detectors,— Si nanocrystal structures,— and two-dimensional transport in laterally 
patterned GaAs quantum wells^. It is remarkable that, despite the microscopic differences 
between these systems, the observed behavior is similar. This naturally raises two questions: 
What are the minimal ingredients needed to observe current branch formation? What is the 
role of heterogeneity in the individual element properties on this behavior? 

In this paper, we address these questions by introducing a nonlinear circuit model consist- 
ing of a voltage-biased series array of N ideal negative differential resistance elements, each 
connected in parallel to a capacitance as shown in Fig. la. We find that heterogeneity in 
the element properties leads to the formation of multiple branches in current-voltage curves 
and a non-uniform distribution of voltages across individual elements. 

The model is defined as follows. Each element of the series array is assumed to have an 
intrinsic I — V curve of the form Ii(vi) with a typical example shown in Fig. lb.— The total 
electrical current is 



for i = 1,2, ...N, where Vj denotes the voltage across the i-th element and is the parallel 
capacitance associated with the i-th element. The total applied voltage to the array is 
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(1) 



V = J2j=i v j- Dividing both sides of Eq. (JTJ by Cj and summing over i allows to express 




where C^J = Y&=i Q" 1 denotes the inverse total equivalent series capacitance. Relabeling 
the summation in Eq. (j2J) as over j and substituting the result back into Eq. (pp), the circuit 
model equations can be written in the following form: 



where = q*q-,- Equation ([3]) has the form of an iV-dimensional dynamical system subject 
to a constraint on the total voltage. The constraint is built into the structure of the model, 
as seen by summing Eq. ([3]) over all % and using the symmetry property = Kji. The 
structure of Eq. ((3]) is equivalent to a limiting case of a standard rate equation model for 
superlattices, when the doping level in each quantum well is very large.- 1 ^ 1 ^ 

In this paper, we focus on the case of linear ramping (with ramp time T) in the total 
applied voltage V, so that V increases linearly between values and NV max and V = 
N^Y^ = Na. The first term in Eq. ([3]) describes the effect of changing total applied voltage, 
while the second term describes a global coupling between individual circuit elements. When 
the ramping rate is large, the effect of the coupling term is expected to be small. On the other 
hand, when the ramping rate is small, the coupling term plays a decisive role in the observed 
behavior. The coupling term displays either positive or negative feedback depending on the 
state of each element. Considering pairs of elements, they attract one another when both 
are on a stable branch of the I — V curve (i.e., Regions I and III in Fig. lb), and repel if 
they are both on the unstable part (Region II). If one element is in a stable region and the 
other is in the unstable region, the elements may either attract or repel depending on their 
relative current values. 

If the coupling term in Eq. fl3]) is neglected, one finds an uncoupled state in which the 
element voltages increase independently. Starting from the initial condition v\ (0) = for 
all i, typical for experimental measurements,-*^ the voltage of each element is given by 
v!> ' — V x Ctot/Ci. Since the C{ are assumed to have a small dispersion, the uncoupled 
state is associated with a field profile that is nearly uniform. (Here, we use the term field to 
refer to the spatial distribution of element voltages.) When the coupling term is nonzero, 
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(3) 
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the system deviates from the uncoupled state and the field profile becomes non-uniform. 
This behavior is usefully characterized by introducing a field inhomogeneity parameter r max 
defined by: 



The quantity r(t) expresses the time-dependent level of field nonuniformity in the system 
for ramp time T, while r max gives the maximal degree of nonuniformity during the entire 
ramping process and associates a single value to the entire process. For the uncoupled state, 
r max tends to zero; however, when the effect of the coupling term is large, r max assumes a 
value of order 1. 

If the total voltage is held constant, the system always relaxes to a state in which the 
currents in the different elements are identical. This corresponds to a stable fixed point of 
the system with V = 0. If the average applied voltage, V/N, falls in the NDR region of 
the single element I — V curve, there are multiple fixed points corresponding to distinct 
arrangements of individual element voltages. For a completely homogeneous system, such 
that all elements are identical, i.e., Cj = C and h{v) = Io(v) for all i, these fixed points 
are degenerate corresponding to the same overall device current.— When the homogeneous 
system is subjected to a ramped voltage, starting from the initial state Vi(t = 0) = for all 
i, the element voltages remain identical throughout the ramp process, i.e., Vi(t) = V(t)/N, 
and the overall I — V curve has a similar shape to that of a single element^. 

When heterogeneity is introduced into the system, current branches emerge in the limit 
of slow ramping as shown in the I — V curve of Fig. 2a. The variance in the individual 
element I — V characteristics is expressed as Uivi) = I (vi)(l + £j), where the e^'s are 
independent and identically distributed (i.i.d.) random numbers, distributed normally with 
mean zero and standard deviation 07, i.e., cr| = 4 Y^Li e l-^ We have also explored the effect 
of capacitance variation and variation due to electrical noise in individual elements, and find 
that the qualitative behavior is the same as that observed when the only variance is in the 
element I — V curves.— 

When the ramp time T is large, the system exhibits well-defined current branches in the 
static I — V curves, shown in Fig. 2a for N = 8. As the total voltage ramps higher, the 
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max 



max{r(t) : < t < T}, 



(4) 



where 
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elements pass very rapidly from Region I to III one by one, and the current exhibits an abrupt 
step down for each such passage.— The parameter r max takes a value of order 1, indicative of 
a non-uniform distribution of element voltages. Figure 2b shows the corresponding contour 
plot in which the current level is plotted in gray scale vs. element number i and total 
voltage V. The clear horizontal bands indicate that, for all parts of the ramping process, 
the individual device currents are identical to one another. 

As the ramp time decreases, the element current levels become somewhat different from 
one another, see Fig. 2c. In this case, more than one element can jump to Region III at the 
same time. The field distribution deviates less from the uniform state, the current branches 
are rounded and smaller in number and amplitude, and the abrupt jumps between current 
branches disappear. A similar rounding of experimental current branches versus ramping 
rate has been reported in weakly-coupled GaAs/AlAs semiconductor superlattices.— The 
r max value also decreases from its large T value. Figure 2d shows the corresponding current 
contour plot in which the horizontal bands are still evident, but interrupted by localized 
dark areas that correspond to the passage of individual or pairs of elements through Region 
II. 

For fast ramping, the element currents are essentially independent of one another, and 
the element voltages pass through the NDR region simultaneously (Fig. 2e). The r max 
parameter approaches zero, implying that the system behavior is very close to the uncoupled 
state described previously. The I — V curve of the full array follows closely that of an 
individual element, Iq(v). The contour plot, Fig. 2f, shows smooth behavior and horizontal 
features are absent. We have also considered ramping from different initial states^ as well 
as more elaborate circuit array models - e.g., including small series inductance and parallel 
capacitance with each nonlinear element^ - and find a qualitatively similar behavior as 
described above. 

Plotting the value of r max versus ramp time T for several different N values in Fig. 3a 
indicates that the transition from uncoupled to fully coupled behavior is a smooth transition 
with an onset that is independent of system size N. For large ramp time, r max approaches 
an asymptotic value that is iV-dependent. As iV — > oo, r max — > 2.85, a value that depends 
only on the shape of the single element function Io(v)r^- For the parameters used here (i.e., 
C = 200 nF and 07 = 0.1), the transition from uncoupled to coupled behavior occurs as the 
ramp time increases from ~ 100 fj,s to ~ 10 ms, a timescale range that is significantly greater 
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than the zero-bias characteristic RC time constant associated with an individual element, 
i.e., t = Cdlo/dv |„=o— 2 ^s. 

To better understand this behavior, we investigate the effect of perturbations about the 
uncoupled state of the form 

Vi = vf 1 ' + 5v{ = v + 5vi, 

where v^(t) = at denotes the uncoupled state solution. Substituting this form into the 
dynamical model, Eq. (J3]), one finds the following system of differential equations valid to 
first order in the e^'s, 



Sv t ~ -i{JoM f>, - *) + I'M J2(5v j - 5 Vi )} = -^^e, - Mfo, ( 6 ) 

where I' Q = dl /dv. Equation (jSj) can be solved explicitly to yield the following expression 
for 5vi\ 

Svi = -—e c a I I ( v ')e c<* dv, (7) 
ua Jo 

with the initial conditions Svi(0) = for all i. 

For large values of ramping rate a, the two exponentials in Eq. (JTj) are approximated by 
unity, and it immediately follows that 

SviCt) = -- t Uat')dt' = -— [ V Uv')dv' = --PM. 
c Jo ca Jo ca 

Substituting this result into the definition of r max , cf. Eqs. (jl]) and ([5|), allows us to write 

n 2 T 2 

■ - 1 -mzx{P{v)/v) 2 



max 

->2 



" tux r 2T/2 

which demonstrates that r max is proportional to T 2 with a coefficient of proportionality that 
is independent of N; this is confirmed in Fig. 3a. 

For somewhat smaller a values, and provided the local maximum of the element I — V 
curve is sufficiently sharp, Eq. (JTj) can be evaluated using a saddle point method to write 
Svi ~ eiI fiL^ e Af c° 0< " , where (5 = —^-4^ \v M > 0. Inserting this result back into the definition 



' (3 ca 

of r(t), cf. Eq. ()5]), gives 

^ , _ KCJjImT 2(I M -I (v)) 

fiCV max v 2 

In the preceding equation, we note that the exponential term takes a maximum value if we 
set v = v m , so that Im — Io( v ) = Im — I m = A/, where v m and I m are the voltage and 



current coordinates of the local minimum in the element I — V characteristic, cf. Fig. lb. 
This allows us to write the maximum value of r over the entire ramp process as 

-e c° . (8) 



Figure 3b plots the values of r max calculated from both Eqs. ([7j) and (jSJ) and compares 
them with the numerical results for iV = 400, demonstrating that there is a range of ramp 
times (around 10 -4 seconds) for which the asymptotic expression, Eq. (jSJ) follows closely the 
transition from uncoupled to coupled behavior. In order to develop an analytic criterion for 
the onset of coupled behavior - and associated field non- uniformity, we define a characteristic 
value of ramp time T t h that sits in this range. Denoting the corresponding characteristic 
value of r max by r th , we have lnr th = 2 In 07 + In (j^r) -lna+ ^f. Solving for T th , we 
can write 

T th = -^pina 7 + K, (9) 
where K = ( \ In r *^ g ^< _ l\ n T th J is slowly varying. Equation (ED shows ex- 



plicitly the dependences of T th on element heterogeneity (i.e., 07) as well as the single 
element NDR behavior. Interestingly, the pre-factor of the In 07 term, calculated to be 
- ^ n:c ~ 1.18 x 10~ 4 s, has the form of an RC time with an effective resistance that depends 
only on the total current drop across the NDR region, AI. In particular, this time scale is 
not sensitive to the shape details of the Io{v) curve in the NDR region. Figure 3c plots T t h 
versus the variance level (with the specific choice r t h = 0.1) and shows good agreement with 
numerical data for N = 40. 

We have introduced a simple model that reveals how negative differential resistance and 
element heterogeneity are key sources for the observation of non-uniform field distributions 
and multiple current branches in the electrical conduction properties of a series array of 
nonlinear circuit elements. Specifically, we have shown how the system approaches a state 
in which the element current levels are fully coupled as the element variance level and 
voltage ramp rate are varied. The numerical and analytical results obtained from this model 
provide insight for understanding similar observed behaviors that are found for a range of 
more complex electronic systems, for example, semiconductor superlattices and quantum 
cascade laser structures.- - - In practical devices it is often desirable to have an electric field 
distribution that is as spatially uniform as possible. The model introduced here points out 
the relevance of two factors for achieving such field distributions. One of these factors is 



the effective level of element heterogeneity in a periodic structure, which, in a device such 
as a superlattice, is determined in part by the fabricated sample quality and the electrical 
noise level. The second factor suggests the possible use of a time-dependent voltage bias 
that possesses a ramping rate through the NDR region of the device that is large enough to 
maintain a nearly uniform field distribution. 
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may be derived for each of the separate current branches. 
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FIG. 1. (a) Nonlinear element array model. Negative differential resistance elements are connected 
in series with parallel capacitance associated to each element, (b) The I — V curve for a typical 
element. Region (II) is the NDR region. 
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FIG. 2. Current-voltage (I — V) curves and current level contour plots for different ramp times, 
all with N = 8, C = 200 nF and 07 = 0.1. (a) Fully coupled case: I — V curve with ramp time 
T = 500 ms. Dashed curve shows the parameter r(t). (b) Current contour plot for T = 500 ms 
(gray scale in units mA). (c) I — V curve for the partially coupled case with T = 2 ms. (d) Current 
contour plot for T = 2 ms. (e) I — V curve for the uncoupled case with T = 0.1 ms. (f) Current 
contour plot for T = 0.1 ms. 
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FIG. 3. (a) The parameter r max versus ramp time T for indicated values of N. Each curve is an 
average of ten different statistical configurations of the Ii(vi) curves with C = 200 nF and aj = 0.1. 
(b) Closed circles give T max vs. T near the characteristic value for N — 400. The solid curve is a 
first-order calculation using Eqs. (jU) and (J7|), while the dashed curve is based on Eq. (jHJ). (c) T t h 
versus 07, comparing analytical prediction with simulation for N = 40. 
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